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1. Introduction 

Finding the density for the first time a Brownian motion reaches a 
moving boundary is an old and classic problem with a wide range of ap- 
plications in both mathematics and applied sciences. The problem may 
be traced back to Louis Bachelier's doctoral thesis [Bachelier (1900)], 
Theorie de la Speculation, and to a paper by the Austrian physicist 
Erwin Schrodinger [Schrodinger (1915)]. 

On the other hand, Robbins and Siegmund (1973) and later Groene- 
boom (1987), Salminen (1988) and Martin-L6f (1998) are among the 
first to explicitly establish the relationship between hitting time prob- 
lems and parabolic partial differential equations. 

The aim of this work is not to provide a new set of solutions to 
the problem of boundary crossing. Our goal is instead conceptual: we 
establish relationships between (a) hitting time problems of Brownian 
motion, (b) 3-dimensional Bessel bridge, (c) Schrodinger's equation 
with linear potential, and (d) heat equation problems with moving 
boundary. We determine conditions under which these problems are 
equivalent. 

The paper is organized as follows: The hitting time problem is in- 
troduced in Section [21 and by use of Girsanov's theorem, we relate it 
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to a functional of a 3-dimensional Bessel bridge. This relationship is 
established in Section [3j 

In Section H] we derive the Feynman-Kac representation that will 
allow us to determine the correspondence between Bessel bridges and 
Schrodinger's equation with linear potential. Next, in Section we 
derive a pair of equivalent problems which in Section O are used to link 
boundary crossing probabilities with the heat equation. We conclude 
in Section [7] with some final remarks. 

2. The problem 

Throughout this paper B stands for a one-dimensional standard 
Brownian motion, and /(•) is a real- valued function denoting a "moving 
boundary". Moreover 

(1) T:=ini{t>0\B t = f(t)} 

is the first time that B reaches the moving boundary /. The main 
motivation of this work is to study the density of T. To be more 
precise, we consider the following problem. 

Problem 2.1. Determine the density ofT, given that, for a > 

(2) f(t):=a+ f f{u)du, /"(£)> 0, and f {f{u)fdu < oo 

Jo Jo 

fort > 0. 
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To analize Problem 12.11 it is convenient to note that using ([T]) and 
setting 



(3) Bt := Bt - [ f\u)du 

Jo 



we obtain the following alternative representation of T: 

T = inf It > 0\B t = a+ / f{u)du\ 
= inf It >0\B t - f'(u)du = a\ 
(4) = inf it > Q\B t = a). 



Since our study will make use of Girsanov's theorem [see for instance: 
Section 3.5, Chapter 3 in Karatzas and Shreve (1991)], let us first 
introduce the heuristics that link these ideas: 

2.1. Girsanov's theorem and hitting time problems. Notice that 
if there exists a measure P under which the process B is a one-dimensional 
standard Brownian motion, then for T as in ([T]) [or (j3j)] and B as in 
©, then 

,-t 



P(T< t) = / h(s,a)ds, t > 0, 
Jo 



where 



(5) h(s,a) := , exp <—-—>, s > 0, a ^ 0. 



II l' 9 

a a 
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This function h is the hitting time density of Brownian motion reaching 
a fixed level a. Alternatively, h is also referred to as the derived source 
solution of the heat equation [see for instance Rosenbloom and Widder 
(1959)]. 

This last statement follows from observing, in (jJJ), that the original 
problem is reduced to that of hitting a 'fixed' boundary a. 

Assuming such measure P exists, how is it related to the original 
measure P? More explicitly, what is the connection (if any) between 
P(T < t) and P(T < £)? 

To answer this question we first need to introduce, the so-called P 
exponential martingales: 

Remark. Exponential martingale and Novikov 's condition. Girsanov's 
theorem states that the process B in ([3]) will be a one-dimensional stan- 
dard Brownian motion under the measure P, as long as the following 
exponential process 



(6) Z t : = exp it f'(u)dB u -If (f'(u)) 2 du\ < t 



< oo 



is a F-martingale. In turn, for this last statement to hold, a sufficient 
condition is Novikov's condition [see for instance, Proposition 3.5.12 in 
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Karatzas and Shreve (1991)] : 



(7) E 



cxp 



(f(u)) 2 du 



< oo VO < t < oo. 



Now, given the P- martingale Z in (Q, together with the P-martingale 



(8)Z t := exp { - f f(u)dB u - i f (f\u)) 2 du } < t < oo, 



where -B is a Brownian motion under P, it follows from Girsanov's the- 
orem that the measures P and P are equivalent and related as follows: 

If (CJ) holds and Z and Z are as in (jH|) and (jHJ), respectively, then the 
pairs (P, B) and (P, S) satisfy the following identities: 



¥[B t eA 



E 



E 



Ztlri, 



(Bt&A) 



Ztl {Bt-S* f'(u)dueA) 



and 



(9) 



F(B t eA) = E 



E 



Zt\B t tA) 



z t \ 



(B t +J*f>(u)dueA) 



This sequence of ideas leads to the following representation: 



Proposition 2.2. Given that T is as in (QJ) and condition |lj) holds, 
the exponential process Z in (G|) «s a P '-martingale. Furthermore, given 
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the hitting time density h in fifty, it follows 



P (T <t) 



(10) 



-/'(*)«-§ /o*(/W*te 



e S°f'{u)B u du 



T = s 



h(s, a)ds, 



fort > 0. 



Proof. The first claim is justified from (j2j) and Novikov's condition, 
equation (J7|). Next, from Girsanov's theorem and the representation 
(J2J), we have the following equality 



P(T < t) = E 



Zt\T<t) 



where Z is as in OH]). 

From the fact that Z is a P-martingale and I(T<t) is J-TAr m easurable, 
it follows from the Optional Sampling Theorem that 



\T<t) = E E(Z t |J- TAt )I 



HT<t) 



(11) 



E 
E 



^tAT I(T<t) 



■^T I(r<t) 



Finally, from (TIT]) together with the integration by parts formula: 



e - £ f{u)dB u -\ ti(f'(u)) 2 du = e -Btf'(t)+f* f"(u)Budu-y*(fi(u))Zdu^ 
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Then the terminal condition Bt = a, and conditioning with respect to 
T under P completes the proof of the proposition. □ 

In view of ([TO]) , our study now focuses on the following expectation: 

Jof"(u)B u du 



E 



T = s 



Thus, our next task is to understand the dynamics of the process B u 
(for < u < s), which is conditioned to hit the fixed level a for the 

first time at time s. 

3. Three dimensional Bessel bridge 

The first idea that comes to mind is that the conditioned process B u 
(for < u < s) is a Brownian bridge, since at t = 0, B = and at 
time s, B s = a. However, this argument is flawed because, in contrast 
to a Brownian bridge, the process B can never be above level a before 
time s. 

Elaborating on the heuristics, let us for the moment write the con- 
ditioned process as: B u . Alternatively, setting 



X u := a — B u , < u < s 



we shall analyze the dynamics of the process X: 

(12) X = a; X u > for < u < s; X s = 0. 
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That is, X is strictly positive (Bessel part), except at time s, when it 
is pinned down at zero (bridge part). It may be formally proved that 
the process just described is in fact the following. 

Definition 3.1. A 3-dimensional Bessel bridge is a stochastic process 
which has the following dynamics: 
1 X, 



(13) dX t :-- 



dt + dW t , X = a>0, 



Xt (s-t) 

< t < s, and W is a one-dimensional standard Brownian motion. 
Furthermore, letting 

(14) k{a,K):=-^exp {-?-], a > 0, K £l 

also known as the source solution of the heat equation [see for instance 
Rosenbloom and Widder (1959)], the transition probability of the pro- 
cess X is: 

G(t,x;r,y) := ¥^{X T edy) 

nK\ V s-tk(s-T,y) 

15 = — ——[k{r-t,y-x)-k{T-t,x + y)] 

x s — t k(s — t,x) 

for < t < r < s, and x, y > 0. 

3.1. Verification of the 3-dimensional Bessel bridge (sketch). In 
order to verify that (fi~2l) and (TT3|) are equivalent, we will make use of 
ideas described in for instance: Imhof (1984) or Chapter XI in Revuz 
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and Yor (2005). First, one derives the finite-dimensional distribution 
of the conditional process (Y t ,0 < t < s\T a = s) with Y t = a — B t . 
[A similar analysis, is carried out for the Brownian bridge in Problems 
5.6.11 and 5.6.13, Karatzas and Shreve (1991).] Later, one may verify 
that Y is a version of the process X, introduced in Definition 13.11 
To this end, let us first recall the following 

Remark. (Brownian motion absorbed at y). Consider the transition 
density of Brownian motion 



G(s,x;t,y) := k(t-s,y-x) 



where k is as in (TTJJ and y < a. It follows — from the reflection 
principle — that the distribution of a Brownian motion, started at x 
at time s and absorbed at level a is given by: 



F s ' x {B t <y,T a >t) = F s '%B t <y)-F^(B t <y,T a <t) 

= F s ' x (B t < y) - ¥ s ' x (B t >2a-y) 
= ¥(B t -B s <y-x) 

-¥(B t -B s >2a-y-x). 
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Or, equivalently: 



¥ s ' x (B t edy,T a >t) = G(s,x;t,y)-G(s,x;t,2a-y) 
(16) = k(t — s,y — x) — k(t — x, 2a — y — x). 



Constructing the finite dimensional distribution ofY. Next, given = 
t < ti < ■ ■ ■ < t n < t n+1 = s, and a positive sequence of values 
Xi, . . . , x n where Xq = a and x n +i = 0, compute: 



P (Y tl e dx u . . . , Y tn e dx n , T a e ds) 

= f{Y ti ed Xi ,l<i<n,T a eds) 

= P (b u ed{a-Xi),l<i<n,T a e ds) 

= P (B u - B ti _ x E d(xi_i - x t ), 1 < % < n, T a G ds 
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where the second equality follows from Y t = a — B t and the third from 
the independence of increments. In turn, this implies that 



F (Y tl e dx u . . . ,Y tn e dx n \T a = s) 

_ P (Y tl e dxi, ...,Y tn e dx n , T a e ds) 
P(T a e ds) 



F(T Xn ed(s-t n ))^ 



F{T a eds) 



IJ \ G (tj-i, a - Xj-ii tj,a- Xj) 



i=i 



•Cr\tj—X) Q> — Xj—l] tjjQ, -\- Xj) 



the last equality follows from equation ( TT6l) (absorbed Brownian motion 
at y), 



'g nh-^-^ 



i=i 
— k(tj — tj-i, Xj + Xj-i) 



Jb y\ O tv \ O ^"Til ** / ' 



a s — t n k(s, 



~\ 11 p(^i tj_i,Xj Xj-x) 



n 



i=i 

— fc(tj — tj_i, Xj + Xj_i) 

Xj_i s — t{ k[s — £j_i, Xj_ij 

x [A; (tj — tj-x,Xj — Xj-x) — k(tj — tj-x, %j + x j-i)\ 
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By inspection one observes that the transition probability of the con- 
ditioned process Y is of the form (TT5|) . For the remaining technical 
details the reader may consult Chapter XI at Revuz and Yor (2005). 

3.2. Hitting problems and 3-dimensional Bessel bridges. From 
Proposition 12.21 and Definition 13.11 we are now able to represent a 
class of first-hitting time problems in terms of a 3-dimensional Bessel 
process, namely: 

Theorem 3.2. Given that h is as in $)), the process X is a 3-dimensional 
Bessel bridge, Definition ^. 1\ Then, the distribution ofT equals 



P(T < t) 



E 



(17) 



<>X P 1 ~ / f"(u)X u du 



1 



x oxp \ -- I (f\u)) 2 du - /'(0)a \ h(s, a)ds 



Proof. Directly verified from Proposition 12.21 and Definition 13.11 □ 

4. Feynman-Kac and 3-dimensional Bessel bridge 

In this section, we will provide a Feynman-Kac representation [see 
for instance: Theorem 4.4.2 in Karatzas and Shreve (1991)] of 



v{t,a) =E 



exp 



f"(u)X u du 



X t = a 
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where X is the 3- dimensional Bessel bridge, equation ( 1T3|) . introduced 
in Definition 13.11 The idea is to apply Ito's rule to the process 

v(y, X y ) expl- f"{u)X u du I , y e [t, s], 

given the following growth and/or boundary condition 



< v(t,a) < 1, < t < s, a> 0. 

We will also make use of the fact that f"(u)X u > for < u < s; that 
the process X is continuous and strictly positive on u G [0, s), and that 
the moments of the running maximum of the process X are known and 
finite. 



Theorem 4.1. Suppose that v(t,a) : [0, s] x IR + — > IR + is continuous, 
of class C 1,2 ([0, s) x R + ),and satisfies the Cauchy problem 

1 d 2 v ( 1 a \ dv 



;i8)- § + /"(*)m; 



f (s, a) 



+ - 



2 9a 2 \ a s — t J da 
1, a G 



ur [°' *) >< 



as well as 



<u(i,a) < 1 V t,ae 



Then v(t, a) admits the stochastic representation: 



v(t, a) = E 4 ' 



t.a 



exp \ - / f"(u)X u du 
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Proof. We proceed as in the proof of Theorem 5.7.6, pp. 366-367, 
Karatzas & Shreve (1991). Applying Ito's rule to the process 



v(y, X y ) expl- f"{u)X u du \ , 
y G [t, s], we obtain, with r n := inf{£ < y < s\X y > n}, 



v(t,a) = E' ,a 
The second term converges to zero as n — > oo, since 



exp < - / f"(u)X u du > I( Tn>s) 
v(r n , X Tn ) exp <^ - / f"{u)X u du \ \ Tn < s ) 



E 



L.a 



v{r n , X Tn ) exp <^ - / f"(u)X u du \ I (Tn < s) 



<E 



t.x 



v(r n , X Tn )I {Tn < s) 



< P r > x (r„ < s) 



F T ' X max X > n 
t<e<s 



E- 



t.x 



< 



max Xfl m 

i<6»<s 



n 



2m 



[see Pitman and Yor (1998) for the moments of the running maximum 
of X] . Finally, the first term converges to 



E'- a 



ox]) <; - / f"(u)X u du 



either by the dominated or by the monotone convergence theorem. □ 
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Our next step is to relate Theorem 14.11 to two alternative, better 
known problems. Namely, that of Schrodinger's equation with linear 
potential. To this end, we first have: 

Proposition 4.2. Solutions to l[T8\) may be of the following form: 



(19) 



v(t, a) 



w(t, a) 
h(s — t, a) 



where 



(20) - w t (t, a) + f"(t)aw(t, a) = -w aa (t, a), [0, s) x 



and h(s, a) is as in |3p. 



Proof. Setting 



and 



u(t, a) = l/h(s — t, a) 



(21) 



v (£, a) := u(t, a)w(t, a), 



we have that 



u t (t, a) 
u a (t,a) 

Uaa(t, a) 



2(s-£) 2 2(s-t) 

1 a 

u{t, a) 



u(t, a) 



a s — t 

1 a 2 

+ 



a 2 (s-t) 2 s-t 



u(t, a) 
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Alternatively, from ( 12 ip it follows 

v t = u t w + uw t v a = u a w + uw a 



u aa w + 2u x w w + uw a 



J aa Uj aa IJJ i ^ j "'x UJ w i u ' UJ aa 



Hence from ffTS]) and fl2T?|) we conclude that 



-u t - -u aa 



as claimed. 



a s — t 



Un 



W 



1 a 

Ua+ [ 

a s — t 



u 



U'n 



D 



The previous Proposition may be generalized by first noting that 



1 x \ hJs — t,x) 



x s — t J h{s — t, x) 



for h as in ([5]). Second, recall that h is a solution to the heat equation 



<H r) '^xx- 



This leads to 



Proposition 4.3. Ifh, v, andw satisfy the following system of partial 
differential equations 



(22) 



i_ . u — )_u 

at * 2 xx 

-v T + k(r, x)v = \v xx + ^f% 
-w T + k(r, x)w = \w xx 
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where a t is only a function of t, then they satisfy the identity: 



(23) 



v(t,x) 



w(t,x) 
h(r, x) ' 



where 



(24) 



r :- 



a u du. 



Proof. Given r is as in (I24p : 



dh Or dh 
~dt ~ ~dt!fr' 



Next, we use (1231 to compute 



w T h — wh T 



w x h — wh x 



h 2 ' x h 2 

(w xx h — wh xx )h 2 — 2hh x (w x h — wh x ) 
v xx = - 4 • 

Finally, substituting in the second equation in f|22|) we have: 



w T h — wh T w 
h? +k h 



1 w TT 1 h TT w 



1 h 2 h 2 
h x (w x h — wh x ) h x 
h? + ~h 



w x h — wh x 
h 2 



which implies 



' h h ~ 2 h 



as claimed. 



□ 
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5. SCHRODINGER'S EQUATION WITH LINEAR POTENTIAL 

Equations ffT9|) and (120]) . together with the growth condition 
(25) < v(t, o) < 1, < t < s, aeR + 

will allow us to introduce a couple of Schrodinger-type problems. 

First Schrodinger problem. From ( 125]) and recalling the definition of h 
in ©: 

II ( 2 

h(t,x) = —? 
we have that 



" (M) = 7!5P exp \"a 



<w{t,a) < h(s-t,a), V0 < t < s, a e R + . 
Furthermore, note that 

lim h(s — t,a) = 

t— >s 

lim h(s — t,a) = S(s). 

The last statement should be interpreted in the following sense: Sup- 
pose A : M + — » M + is a C 1 function, then 

lim / h(s - t, x)X(t)dt = X(s), < t < s < oo. 

This leads to our first representation: 
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Remark 5.1. Let v be as in Theorem \4-l\ and h as in |3]). Then the 
following relationship holds 



w(t,a) r , m4 

v(t,a)= K \ J te[0,s]xR- ] 

h(s — t, a) 



Furthermore w is a solution to the Schrodinger equation with linear 
potential 



(26) -fto(t, a ) + f{t)aw{t,a) 



1 d 2 w 

2 0^*'^' te ^' aG 



satisfying the following boundary conditions 



limw(t,a) = 

t— >s 



lim w(t, a) = S(s), 

a->0 



as well as the compatibility condition 

t-*s h(s — t, a) 
For our second problem, we will make use of the following: 



-Wv 
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Proposition 5.2. If h, k, v, and w satisfy the following system of 
partial differential equations 

_J_ h — lh 

— h — -h 

-v T + a(r, x)v = \v xx + !fv x 

+ [a(r,x) + (£-£)£]«,= \ Wxx + «;«(%■- ^) 
where at is only a function of t, then they satisfy the identity: 

w(t,x) 

v ( r > x ) = T? \> 

h{T, x) 

where 



ft 
(27) r :— I a u du. 



Proof. We proceed as in the proof of Proposition 14.31 □ 

Second Schr 6 dinger problem. Given the fundamental or source solution 
of the heat equation k, equation (TT4|) . we may relate v to another 
Schrodinger-type problem, by setting 



. . u(t,a) 

v(t,a) ~ 



k(s — t, a) 
and making use of Proposition 15.21 
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Remark 5.3. Let v be as in Theorem \4-l\ and k as in equation HW , 
then the following relationship holds: 



v(t, a) 



u(t, a) 



k(s — t,a) 



t e [0, s) x R^ 



Furthermore, it follows from Proposition \5.2\ that u is a solution to the 
Schrodinger equation with linear potential: 



(2* 



du 
"dt 



{t,a) + 



1 



(s-t) 



+ f"(t)a 



u(t, a) 



1 u 1 du 

-—(t,a)+- ir (t,a), te[0,s], aeR + 

2 aa A a oa 



satisfying the following boundary conditions: 



limu(t, a) = 5(a) 



limw(£,a) = 0, 

a— >0 



as well as the compatibility condition 



u{t,a) 

hm-- 

t-*8 k(s — t, a) 



1. 



The boundary conditions of this second problem follow, once more, 
by using the growth condition in Theorem 14.11 i.e. 



0<u{t,a)<k{s-t,a) Wte[0,s],ae 
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We conclude this section by making a few remarks. First, note that 
the problem described in Remark 15.31 is a Cauchy problem, the one 
appearing in Remark 15. II is not. 

We have introduced the consistency conditions in Remarks 15.11 and 
15.31 in order to make them compatible with Theorem 14.11 Boundary 
conditions for the density of T will be derived in the following section. 

The next observation is that for the equations ( 126]) and (128]) one 
may construct particular solutions. There is in fact one solution of 
( 12"B"]) which may be used to find a new set of solutions of f ]2"5]) . To do 
so, one may use a technique introduced in Bluman and Shtelen (1996). 
However, solutions to the problems described in Remarks 15.11 and 15.31 
for arbitrary bounds / — which do not involve integral equations — have 
yet to be found. 

This last remark will be clarified in the following section. The idea 
stems from the following fact (yet to be shown): There is a particular 
set of solutions of either ( 126]) and/or ( 128]) which transform our original 
problem — that is, finding the density of hitting times — into that of 
finding particular solutions to the heat equation killed at a moving 
boundary! Unfortunately, in practice this is not very useful to us, since 
the problem leads back to Volterra integral equations as described in 
De Lillo and Fokas (2007). 
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A final curious remark is that equation ( T2T)j) may be transformed into 
a non-homogeneous Burgers' equation. The homogeneity is only time 
dependent and it corresponds precisely to /": 



(29) Kt(t, a) + n(t, a) ■ K a (t, a) - -K aa (t, a) = f"(t). 



Alternatively, k and w — equation (|2T)|) — introduced in Remark 15.11 are 
related through: 



d 
K(t,a) = -—logw(t,a) 



From Proposition 14.31 it follows: 

Proposition 5.4. // h, v, w and u satisfy the following system of 
partial differential equations 

— J- . h ■ - ih 

at ' 2 xx 

-v T + k(r, x)v = \v xx + !fv x 
—w T + k(r, x)w = \w xx 



-U T ~r U • U x r,U xx I rC x 



where a t is only a function of t, then they satisfy the identities: 



V[T,X) 



w{t,x) 
h(r, x) 



W; 



u 



w 



w = e-f udx } 



HITTING TIMES, BRIDGES, AND SCHRODINGER'S EQUATION 25 

where 



ft 
t := / a u du. 



Of course, at first sight, this is not so surprising since this relationship 
is well known. The possibly "interesting" part lies in the fact that so 
many probability-related papers have been written regarding Burgers' 
equation (129]) . see for instance Barndorff-Nielsen and Leonenko (2005) 
and references therein. 



6. Particular solutions to Schrodinger's equation with 
linear potential (heat equation) 

Particular solutions to Schrodinger's equation with linear potential 
have been obtained, for instance, by Feng (2001) an references therein. 

The strategy employed in Feng (2001) is not only interesting but it 
is also useful. This follows since it ultimately leads to heat equation 
problems. In this section, we propose yet another equivalent approach 
using Fourier transforms: 
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Theorem 6.1. Solutions to / fUZTj) are given by 



(30) 



w(t,a) = e H s (/») 2a!u + a /'«_ / u(y)e'^ y2{s ~ t)+iy(a+ ^ f ' {u)du) dy 

2tt J_ 00 



e y t a (f'(u)fdu+af(t) 



u I s — t, a + / f'{u)du J 



where u is a solution to the heat equation and IT is an arbitrary func- 
tion, as long as the integral is well defined. 

Proof. Let 



oo 



(31) w(t,X) := / e~ lXa w(t,a)da. 



Applying the Fourier transform to (|20|) we have 



-w t {t, A) + if"(t)w x (t, A) + ^X 2 w(t, A) = i := V^l- 



Next, set y = A + if it) and w(t, A) = u>(£, y), that is: 

W t = Wf + if"(t)Wy, W X = Wy 

which after substitution in ( 13TJ) leads to 



(32) - w t (t, y) + ^(y- tf(t)) 2 w(t, y) = 0. 
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Consequently, setting Ii(y) = w(0,y) it follows that 



w{t,y) = U(y)expl-- (y-if\u)fdu 

(33) = n(y)exp|-iy 2 ( S -t)+iy j° f{u)du+ 1 ~ J\f{u)) 2 du 



[Note that equation ( 1331) will be a solution of ( 1321) for arbitrary IT.] 
Alternatively, this implies that 

w (t f a) = ± / Yi(y) e -b 2 ^-^yS:f'{u)du+\s:{f { u)f d u el xa dX 



2vr J _ 00 

e ^/ t S (/'(«)) 2 ^+a/'(t)J_ /"°° n ( 2/ ) e -^ 2 (^-<)+i?/(a+/ t s /'Hrf«) rf2/ 



as claimed. D 

This leads to our already anticipated claim, that boundary crossing 
probabilities, are equivalent to a set of heat equation problems. 

Theorem 6.2. The density ip of the first passage time T defined in (TJP 
is related to a solution of the heat equation u as follows 



<p(s, a) = u) ( s, a + / f'{u)du 
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and it is bounded by 



h(s, a)e" a/ ' (0) "5 So U'{u)?du-j° /" ( u )e« (x u )du 



< oo [ s, a + / f'(u)du 



</*(*, a)e^°H Jo (/») 2 <^ 



where h is as in |3]]. 



Proof. It follows from Jensen's inequality, 



oxp <{ - J f"{u)W' a (X u )du } < E*' 



t.a 



<>X P 1 - / f"{u)X u du 



< 1, 



equation (1201), and (1301). 



□ 



As a corollary we obtain the necessary boundary conditions for den- 



sity <p: 



Corollary 6.3. From the previous Theorem, the following boundary 
conditions hold: 



(34) 



lim u I s, a + / f'(u)du 

a^O \ Jo 

lima; ( s, a + / f'(u)du 
s ^° V Jo 



6(0), 
0. 
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and 



1 



V27TS 3 



-.e 



S°U'{u)?du-2^l J S f»{u)V5=Zy/Zdu 



< lim uj a I s, a + / f'(u)du ) 
a ^° V Jo J 

< x e -§ r '(/») a du 

a/27TS 3 



or given t/ie definition of Fractional Integral 



1 /' x 
■^(s) := ^7-^ / (a; - y) a ' 1 g(y)dy, x > 0, a e 
1 (a) Jo 



we /iave 



1 e -i ro(/'w) 2 ^-s- l/2 J 3/2 /"(s)y2 

\/2' n " e3 



OTS"' 



< lim uj a I s, a + / f{u)du J 



< _J e -!/o(/») 2 ^ 

V27TS 3 



That is, Theorem 16.21 together with Corollary 16.31 uniquely charac- 



terize the density of T. Alternatively, from Theorem 16.21 we have 
obtained — for free — lower bounds for the problems described in Re- 
marks 15.11 and 15.31 
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7. Concluding Remarks 

As was already mentioned in the Introduction, the aim of this pa- 
per is not to provide a new set of solutions to the problem of bound- 
ary crossing. Our goal is instead conceptual, in that we first make 
use of Girsanov's theorem to connect hitting time problems to Bessel 
bridges. We do so by use of the Optional Sampling Theorem. Next, 
we relate exponential functionals of 3 dimensional Bessel bridges to 
Schrodinger's equation with linear potential. This is done by uniquely 
relating both problems through a Feynman-Kac representation. We 
conclude by characterizing the density of T in terms of a heat-equation 
problem and/or two alternate Schrodinger-equation problems. 

Girsanov's link is useful since it suggests how to extend this idea to 
higher dimensions or more general diffusions. Alternatively, the con- 
struction of Bessel bridges points to the fact that the concept can be 
extended. From the derivation of our Feynman-Kac representation we 
have learned that in some particular cases, sufficient conditions can be 
relaxed. Finally, Schrodinger's equation with linear potential may be 
the way in which boundary crossing probabilities of some "higher or- 
der" processes can be described. For example, instead of having a para- 
bolic equation — which describes the transition probability of Brownian 
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motion — as in (|26p . it could be that a process with transition proba- 
bility 

f)0 r) n f) 

— (£,x) = a— (£,x), [0,oo)xR, n > 2, 

has a first passage time probability which might be described by a 
Schrodinger type-problem: 

Ow n w 

— (t, x) + xf"(t)w(t, x) = a-^{t, x), [0, oo) x R. 

In particular note that the techniques used in the proof of Theorem 16.21 
still hold. 
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